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Abstract 

Self -seeding schemes, consisting of two undulators with a monochromator in be- 
tween, aim to reduce the bandwidth of SASE X-ray FELs. We recently proposed 
to use a new method of monochromatization exploiting a single crystal in Bragg- 
transmission geometry for self-seeding in the hard X-ray range. The obvious and 
technically possible extension is to use such kind of monochromator setup with 
two -or more- crystals arranged in a series to spectrally filter the SASE radiation at 
two -or more- closely-spaced wavelengths within the FEL gain band. This allows 
for the production of doublet -or multiplet- spectral lines. Applications exist over a 
broad range of hard X-ray wavelengths involving any process where there is a large 
change in cross section over a narrow wavelength range, as in multiple wavelength 
anomalous diffraction techniques (MAD). In this paper we consider the simultane- 
ous operation of the LCLS hard X-ray FEL at two closely spaced wavelengths. We 
present simulation results for the LCLS baseline, and we show that this method can 
produce fully coherent radiation shared between two longitudinal modes. Mode 
spacing can be easily tuned within the FEL gain band, i.e. within 10 eV. An inter- 
esting aspect of the proposed scheme is a way of modulating the electron bunch at 
optical frequencies without a seed quantum laser. In fact, the XFEL output intensity 
contains an oscillating "mode-beat" component whose frequency is related to the 
frequency difference between the pair of longitudinal modes considered. Thus, at 
saturation one obtains FEL-induced modulations of energy loss and energy spread 
in the electron bunch at optical frequency. These modulations can be converted into 
density modulation at the same optical frequency with the help of a weak chicane 
installed behind the baseline undulator. Powerful coherent radiation can then be 
generated with the help of an optical transition radiation (OTR) station, which have 
important applications. In this paper we briefly consider how the doublet struc- 
ture of the XFEL generation spectra can be monitored by an optical spectrometer. 
Furthermore, the OTR coherent radiation pulse is naturally synchronized with the 
X-ray pulses, and can be used for timing the XFEL to high power conventional 
lasers with femtosecond accuracy for pump-probe applications. 
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1 Introduction 



LCLS began routine user operation with photon energy up to 10 keV |Ht]. 
This success motivated the planning of a significant upgrade over the next 
several yearsjH]. Plans for the near- term upgrade include self-seeding for 
hard X-rays ||2|]. 

In li-lll, self-seeding options for the LCLS baseline were investigated using 
a scheme relying on a single-crystal monochromator in Bragg-transmission 
geometry The Bragg crystal reflects a narrow band of X-rays, resulting in 
a ringing within the passband in the forward direction, which can be used 
to seed the second undulator. The chicane creates a transverse offset for the 
electrons, washes out previous electron-beam microbunching, and provides 
a tunable delay of the electron-bunch with respect to the radiation, so that 
the electron beam only interacts with the ringing tail of the X-ray pulse 
01-Bl. 

In this paper we extend our previous investigations to the case of a self- 
seeding scheme with a similar kind of monochromator, consisting of two 
or more crystals in the Bragg geometry arranged in a series. The impinging 
SASE radiation is therefore spectrally filtered at two or more closely spaced 
wavelengths within the FEL gain band. This allows for generation of doublet 
spectral lines. 

Simultaneous operation of a high-gain FEL at two separate wavelengths 
using narrow band seed lasers was previously studied in jjj. Applications 
of such mode of operation involve any process where there is a large change 
in cross section over narrow wavelength range, as in multiple wavelength 
anomalous diffraction technique [9(]. Here we present simulation results for 
the LCLS baseline and we show that this method can produce a doublet 
structure in the LCLS hard X-ray FEL spectrum. 

The XFEL-output intensity contains an oscillating "mode-beat" component 
whose frequency is related to the frequency difference between the pair of 
longitudinal modes considered. With the help of a weak chicane installed 
behind the baseline undulator, see rfioll , a powerful, coherent pulse of radi- 
ation can then be generated with the help of an optical transition radiation 
(OTR) station, which has important applications. In this paper we briefly 
consider how the generation of doublet structures in the XFEL spectrum 
can be monitored by an optical spectrometer. Furthermore, the OTR coher- 
ent radiation pulse is naturally synchronized with the X-ray pulses, and 
can be used for timing the XFEL to high power conventional lasers with 
femtosecond accuracy for pump-probe applications. 

1 Corresponding Author. E-mail address: gianluca.geloni@xfel.eu 
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self -seeding scheme with wake monochromator for the LCLS baseline undulator 



two-color operation of LCLS baseline at closely spaced wavelengths 
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Fig. 1. Design of the LCLS baseline undulator system for generation of doublet 
spectral lines. 

2 Operation of the LCLS hard X-ray FEL at closely spaced wavelengths 
by use of a self -seeding scheme with a wake monochromator 



In the present Section we introduce our technique. An overview is sketched 
in Fig. 1. The reader may recognize similarities with the previously proposed 
schemes 101-10], except for the presence of two crystals, instead of one, 
located in the transverse offset provided by the magnetic chicane. As before, 
the first undulator in Fig. 1 operates in the linear high-gain regime starting 
from the shot-noise in the electron beam. 

After the first undulator, the output SASE radiation passes through the 
monochromator, consisting of a series of two (or more) crystals in the Bragg 
transmission geometry. The monochromator setup is sketched in Fig. 2. Ac- 
cording to our monochromator principle, the SASE pulse coming from the 
first undulator impinges on a set of two or more crystals in Bragg diffraction 
geometry. Then, the crystals operate as bandstop filters for the transmitted 
X-ray SASE radiation pulse, as shown in Fig. 3. The incident angles of the 
first and of the second crystal are different, so that the Bragg condition 
is met for different frequencies within the SASE spectrum. The temporal 
waveform of the transmitted radiation pulse shows a long tail, whose du- 
ration is inversely proportional to the bandwidth of the absorption line in 
the transmittance spectrum. This tail can be used for seeding the electron 
bunch after the chicane. It should be noted here that the ringing tail does 
not depend on the distance between the crystals, and that such distance has 
no influence on the output characteristics of the radiation. At variance with 
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operation of wake monochromator at two closely spaced wavelengths 
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Fig. 2. Forward diffraction in a series of crystals in Bragg geometry. The incident 
angles on the first and on the second crystal are different. Due to multiple scattering, 
the transmittance spectrum of the crystals shows doublet absorption lines. The 
temporal waveform of the transmitted radiation pulse is characterized by a long 
tail. The radiation power within such tail is shared between two harmonic waves 
with slightly different wavenumbers. 

previously considered schemes relying on a single crystal, here two band- 
stop filters are present at different frequencies. As a result, both frequency 
components are present in the seeding signal, and will be amplified in the 
output undulator. This allows for the production of doublet spectral lines, 
Fig. 4. The output signal from our setup is thus fully coherent, but shared 
between two separate longitudinal modes. The relative alignment tolerance 
of the crystal tilting angles should now be in the order of ten microrad, in 
order to allow for a stable frequency difference within the Darwin width. 

While the radiation is sent through the crystals, the electron beam passes 
through a magnetic chicane, which accomplishes three tasks: it creates an 
offset for the crystals installation, it removes the electron microbunching 
produced in the first undulator, and it acts as a delay line for the imple- 
mentation of a temporal windowing process. In this process, the magnetic 
chicane shifts the electron bunch on top of the monochromatic tail created 
by the bandstop filter thus temporally selecting a part of it. By this, the elec- 
tron bunch is seeded with a radiation pulse characterized by a bandwidth 
much narrower than the natural FEL bandwidth. For the hard X-ray wave- 
length range, a small dispersive strength R 56 in the order of a few microns 
is sufficient to remove the microbunching generated in the first undulator 
part. As a result, the choice of the strength of the magnetic chicane only 
depends on the delay that one wants to introduce between electron bunch 
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Fig. 3. Two crystals bandstop filter as monochromator for the self-seeding scheme. 
The two crystals in Bragg geometry are arranged in a series, and act as bandstop 
filters for the transmitted X-ray SASE radiation pulse. Each crystal reflects a narrow 
band of X-rays, resulting in a ringing in the pass band that can be used to seed the 
second undulator at two closely spaced wavelengths. 

generation of doublet spectral line of the LCLS hard X-ray FEL at low-charge mode of operation 



first mode second mode 




8 keV photon energy 



Fig. 4. The doublet structure of the X-ray FEL generation spectra. 

and radiation. The optimal value amounts to R 56 12 fim for the low-charge 
mode of operation. Such dispersion strength is small enough to be gener- 
ated by a short (4 m-long) magnetic chicane to be installed in place of a 
single undulator module. Such chicane is, at the same time, strong enough 
to create a sufficiently large transverse offset for installing the crystals. 
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Fig. 5. Sample of out put data as a function of the position inside the bunch. (Upper 
left) Output power after 14 undulator cells. The dotted line represents the electron 
bunch profile. (Upper right) Average energy loss. (Bottom left) rms energy spread. 
(Bottom right) Current profile after a weak chicane installed after the LCLS baseline 
undulator, R56 = 0.08 mm. 



Since two bandstop filters are present at slightly different frequencies, the 
seed signal contains a "mode-beat" component whose frequency is related 
to the difference between central frequencies of the bandstop filters. As a 
result, the output from our setup exhibits a remarkable feature: namely we 
have modulation of the output radiation power and of the energy loss and 
energy spread of the electron beam at a frequency related with the difference 
between the central frequencies of the bandstop filters. This is exemplified 
in Fig. 5, upper right, upper left and bottom left plots, obtained with the 
code Genesis [HQ . In the example case shown in Fig. 5, this modulation is 
around 200 nm. 



One may take advantage of such energy modulation by transforming it into 
a density modulation, which can be used to produce powerful pulses of 
coherent radiation in the visible range. This can be done with the help of 
an extra magnetic chicane installed after the XFEL undulator at the LCLS 
baseline, as depicted in Fig. 6. A dispersive strength of 50 - 100 microns 
is enough to this end. The output particle file from Genesis is used as 
input to the code Elegant to study the electron beam dynamics through 
the magnetic chicane. Tracking results are shown in Fig. 5 bottom right. 
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simplest setup for producing laser lite optical pulse from electron 
bunch when 1XUS hard X-ray FEL generates doublet spectral line 
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Fig. 6. Installation of a small dispersive element (magnetic chicane) and of an OTR 
station after the LCLS baseline undulator will allow to produce laser-like optical 
pulses from electron bunches when the LCLS hard X-ray FEL generates doublet 
spectral lines. 

As a simple example of radiator one may consider an OTR station, as 
sketched in Fig. 6 and Fig. 7. It is technically possible to avoid interference 
between the OTR pulse production and the X-ray beam delivery. In fact, 
the vacuum chamber downstream of the main undulator has an effective 
aperture of 25 mm, which is the diameter of the BPM bore. Such aperture is 
sufficiently large to alter the electron beam trajectory (with the help of cor- 
rectors) and introduce an offset in the order of a centimeter, thus separating 
the electron beam and the X-ray. Since we deal with optical wavelengths, 
such electron orbit perturbation will not perturb the density modulation in 
electron bunch. 

Several applications can be considered. With reference to Fig. 7 the presence 
of coherent pulses of OTR can be used, with the help of an optical spec- 
trometer, to monitor the generation of the doublet spectral lines. Moreover, 
the OTR pulses are naturally synchronized with the electron bunches, and 
therefore with the main X-ray FEL pulses. After filtering, these pulses can 
be cross-correlated with a harmonic signal from a TiSa pump laser in order 
to provide a timing of pump and probe pulses with an accuracy within a 
few femtoseconds. 
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Fig. 7. Block diagram view of possible coherent OTR pulse applications. 
Table 1 

Parameters for the low-charge mode of operation at LCLS used in this paper. 



Units 



Undulator period 


mm 


30 


K parameter (rms) 




2.466 


Wavelength 


nm 


0.15 


Energy 


GeV 


13.6 


Charge 


nC 


0.02 


Bunch length (rms) 


fim 


1 


Normalized emittance 


mm mrad 


0.4 


Energy spread 


MeV 


1.5 



3 FEL simulations 



Following the previous introduction to our technique, in this Section we 
report on a feasibility study performed with the help of the FEL code GEN- 
ESIS 1.3 [11] running on a parallel machine. In this Section we will present 
the feasibility study for the short-pulse mode of operation of the LCLS. Pa- 
rameters used in the simulations for the short pulse mode of operation are 
presented in Table 1. We present a statistical analysis consisting of 100 runs. 

Our starting point is the simulation of the SASE radiation characteristics at 
the double crystal. The incident power distribution after 11 section is shown 
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Fig. 8. Power distribution after the SASE undulator (11 cells). Grey lines refer 
to single shot realizations, the black line refers to the average over a hundred 
realizations. 

in Fig. 8. The bandstop-effect of the two diamond crystals (C(400) reflection, 
100 /.tm-thick crystal) is evident by inspection of Fig. 9. The effect in the time 
domain is shown in Fig. 10. As shown before in [3]-|4], the presence of a 
single bandstop-filter is related with a long, monochromatic tail following 
the main SASE pulse. Such structure has a duration comparable with the 
inverse of the filter bandwidth. Moreover, since there are now two frequency 
filters centered at different central frequencies, a beating between the two 
waves at slightly different frequencies takes place, yielding a seeding power 
distribution modulated at optical frequencies, as is evident from the right 
part of Fig. 10, which is an enlargement of the average (black) line in the left 
part of the same figure. 

Photons and electrons are recombined after the crystals, so that the electron 
bunch temporally selects a part of the long tail in Fig. 10, and is seeded 
by it. The seeded electron beam goes through the output undulator, Fig. 1, 
which is sufficiently long (14 cells) to reach saturation. The output power 
and spectrum of the setup is shown in Fig. 11. The effect of the double 
crystal is clear. The output power consists of a 20 GW pulse of coherent 
radiation, divided in two separate longitudinal modes, yielding doublet 
spectral lines on the right part of Fig. 11. By inspection of Fig. 11 one can 
see that the power is not evenly distributed between the two modes. In 
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Fig. 9. Spectrum after the diamond crystals, 100 jUm-thick, C(400) reflection. The 
bandstop effect is clearly visible. Grey lines refer to single shot realizations, the 
black line refers to the average over a hundred realizations. 

particular, the low frequency mode is slightly dominant. The way the two 
modes compare with respect to each other depends in a complicated way 
on the position, in frequency, of the band-stop filters with respect to the 
impinging SASE spectrum, i.e. on the angular tilt of the crystals. Fig. 12 
shows several simulation cases demonstrating different possible situations. 

It should be noted that the phase of the sinusoidal envelope of the seed signal 
(see Fig. 10 right) is randomly distributed from to 2n on a shot-to-shotbasis. 
As a result, the output power in Fig. 11, averaged over 100 shots, does not 
exhibit any "mode-beat" component. Small residual oscillations, see Fig .11 
left, are only due to the fact that the averaging process is performed on a 
finite number of shots (one hundred) only. 

Additionally, it is interesting to remark that the average output power at 
saturation (15 GW) is about the same that we had in the previous case 
where a single spectral line was generated. In other words here we only 
have a redistribution of the total number of photons between two modes, 
but the total number of photons remains about constant. 

Finally, we note how the spectra in Fig. 12 right, exhibit small sidebands on 
the left, and on the right of the doublet. This effect is ascribed to nonlinear 
interaction of the two seed waves. This leads to generation of a discrete 
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Fig. 10. (left) Power distribution after the diamond crystals. The monochromatic 
tail due to the transmission through the bandstop filters is now evident on the left 
of the figure. Grey lines refer to single shot realizations, the black line refers to the 
average over a hundred realizations, (right) An enlargement of a single shot seed 
power distribution presented in the left part of the plot. The seed signal consists 
of two harmonic waves with slightly different wavenumbers. The form of this 
type of signal is interesting, and consists of a rapidly oscillating (at X-ray angular 
frequency) carrier wave multiplied by a simple sinusoidal envelope function at 
optical angular frequency. The beat of the sinusoidal envelope (see the maximum 
at 6 fira) is related to specific features of the transmittance profile of the band-stop 
filter. 



Fig. 11. (left) Power distribution of the X-ray radiation pulse at saturation, (right) 
Spectrum of the X-ray radiation pulse. Grey lines refer to single shot realizations, 
the black line refers to the average over a hundred realizations. 

spectrum of beat waves extending beyond the FEL gain band. In particu- 
lar, in the spectral window shown in the plots, we can distinguished two 
nonlinear resonances at 2a)\ - coi and 2cu 2 - co\. Generation of such wave 
beat resonance has been demonstrated and studied in the case of LCLS 
parameters in |@]. 
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Fig. 12. Effect of different position of modes within the FEL gain band at fixed mode 
spacing, (left) Position of the frequency filters relative to the X-ray spectrum, (right) 
Spectrum of the X-ray radiation pulse. Grey lines refer to single shot realizations, 
the black line refers to the average over a hundred realizations. 

4 Coherent OTR from an optically modulated electron bunch 



At the exit of the output undulator, the electron bunch exhibits substantial 
energy modulation at wavelengths in the range of the spacing between the 
doublet lines. This behavior, confirmed by simulations, see Fig. 5 bottom 
right is to be expected on the basis of qualitative reasoning as well. In fact, in 
the output undulator microbunching develops at both seeding wavelengths, 
leading to a beating with frequency related to the difference frequency. This 
beating is evident from the analysis of the output power from a given single- 
shot simulation, e.g. the one in Fig. 5, upper right. It follows that electrons 
lose energy following such beating pattern, thus leading to the development 
of a strong energy modulation. 



According to the scheme described in Section 2, Fig. 6, a weak dispersive 
element transforms such energy modulation in a substantial density mod- 
ulation. A longitudinal dispersion of several tens microns is enough to 
optimize the modulation level. Starting with the output from Genesis, we 
performed a simulation of the bunching in the chicane with the help of the 
3D tracking code Elegant [12]. The result is shown in Fig. 5, bottom right 
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plot, where we consider a case study with the beam current profile I(t) after 
a chicane with R 56 = 0.08 mm. Such chicane can be as short as a few meters. 
In this Section we estimate the number of optical photons emitted by the 
modulated electron bunch as it travels between the last bend of the chi- 
cane and the OTR screen, Fig. 6. We base our investigation on the analysis 



presented in [13Q 



Let us define the slowly varying envelope of the field in the space-frequency 

domain as E = Eexp [—icoz/c]. We will refer to this quantity simply as 
■± 

"the field". Here E(co,r,z) is the Fourier transform of the electric field 
_> -± 

E(t, r, z) in the space-time domain, according to the convention: E(co, r, z) = 
E(t, r, z) exp[icot]dt. Note that r indicates the transverse position vector. 

The energy radiated by the bunch per unit frequency interval per unit 



surface can be calculated as in 114J], and expressed in terms of number of 
photons, yielding: 



dN ph _ c 
dcodS An 2 hoj 



E(z,r) 



(1) 



In order to obtain the field produced by the electron bunch at the posi- 

tion of the OTR screen, E(r), a microscopic approach can be used where 
the single-electron field is averaged over the six-dimensional phase-space 
distribution of electrons. In this (Lagrangian) approach particles are labeled 
with a given index, and the motion of individual charges is tracked through 
space. One follows the evolution of each particle as a function of energy 
deviation by, angular direction if, position I and arrival time t at a given 
longitudinal reference-position. Knowing the evolution of each particle, in- 
dividual contributions to the field are separately calculated and summed 
up. Due to the high-quality electron beams produced at XFELs (highly colli- 
mated and nearly monochromatic) we have the simplest possible situation. 
Namely, when performing OTR calculations from an optically modulated 
electron bunch we can neglect both angular and energy distribution, and 
use a model of a cold electron bunch with given longitudinal, f T (t), and 

transverse, //(/), charge density distributions. Note that in general we have 
no factorization of the charge density distribution into longitudinal f T (t) and 

transverse fi(l) factors. However, here we will be interested in an estimate 
of the number of available photons and, with some accuracy, we can use 
a model with separable charge density distribution function. In fact, this 
assumption is not related to fundamental principles, and will only lead to a 
different numerical factor in the estimation of the number of photons. 

Let us first consider the transverse charge density distribution //(/). In the 
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Fig. 13. 1(t)—Io(t) as a function of time. This plot shows that part of the electron beam 
current responsible for electron beam density modulation at optical wavelengths. 



following we will simply assume a Gaussian shape: 



//(0 = 



exp 



2ol 



(2) 



Let us now turn to the longitudinal charge density distribution f T (t). In our 
study case, the electron beam exhibits a relatively broadband modulation 
around the nominal wavelength of 200 nm. It is therefore convenient to 
write f T (t) as 



Mt) = U(t) + fAt) ■ 0) 

Here / x o is related to the unmodulated electron beam current by 



/to(0 = 



1 jog) _ 1 



exp 



24 



(4) 



where L 



3 kA, and 07 — 3 fs is the electron bunch rms duration, 



which we assumed to be Gaussian throughout all this paper. Using a similar 
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definition, the modulation part of the charge density distribution, / T i(f) can 
be found as 



M(t) = — , (5) 

y2na T I max 

and is explicitly shown in Fig. 5 for the numerical case under investigation. 

Having said this, the field distribution for the electron bunch at the OTR 
screen in the space-frequency domain is essentially a convolution in the 
space domain of the temporal Fourier transform of the charge density dis- 
tribution and the temporal Fourier transform of the single-electron field, 
yielding 



E(r) = N e f T (co) J dr'E {r')0- f) , (6) 

where f T (co) is the Fourier transform of the temporal charge density distri- 
bution f T (t), and E is the electric field from a single electron. 

Since we need to calculate the total number of OTR photons available, the 
longitudinal position of the observation plane, z, does not matter, and it is 
easier to perform calculations in the far zone for calculation purposes only. 
From Eq. (6), from the definition of spatial Fourier transform, and Eq. (7) 
we obtain 



dN ph dN® 
— — = N — — 

dcodS e dcodS 



Q(co,Z) , (7) 



where is the number of electrons in the bunch, (1) refers to a single 
electron, and 



g (co,Z)=Mco)Mg), (8) 

fi{£) being the Fourier transform of transverse charge density distribution, 

with £ the angular variable conjugated to r. Let us study each of the factors 
in Eq. (7) separately. 

The spatial Fourier transform of //, i.e. ft may simply be written as 
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Fig. 14. Squared modulus of the Fourier transform of f T \, |/ T il 2 , as a function of the 
frequency. 



//(£) = exp 



a^a 2 
2c 2 



(9) 



which is dimensionless. / T can be calculated as the sum of two parts: 



AM =/toM+/tiM, 

where 



(10) 



/ T oM = exp 



2 2 
OjCO 



(11) 



and / T i is calculated numerically, |/ T i| 2 being shown in Fig. 14. 



Finally, the photon density distribution from a single electron can be cal- 
culated from the knowledge of the field angular distribution relative to a 
single electron according to: 



^ c 

dcodS An 2 fyi) 



E (z,£) 



(12) 
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Following |13|], one obtains the far-zone expression for the single-particle 
field 



E(z,l) = 



cz(y 2 l 2 + 1) 



exp 



2c 



(13) 



As said before, Eq. (13) is the spatial Fourier transform of the field at the 
OTR screen: 



c l y r \cy ) 



(14) 



where K\ is the first order modified Bessel function of the second kind. Eq. 
(14) is valid, strictly speaking, only in the region r « yA on the OTR screen. 
However, as noted in [13j], corrections to Eq. (14) enter only logarithmically 
in the calculation of the number of photons outside such region. Therefore, 
as concerns photon number estimations, we may still use Eq. (14), yielding 
Eq. (13) which in its turn, substituted in Eq. (12), gives 



_ p ± = -2L l± d5) 

dwdS n 2 wz 2 (y 2 e + l) 2 ' K } 

with a = e 2 /(hc) = 1/137 the fine structure constant. 

Substitution of Eq. (15) and Eq. (8) in Eq. (7) should finally be followed by 
integration over a proper bandwidth and angular range to yield the number 
of photon available. 

Application, with some accuracy, of the adiabatic approximation, meaning 
that the bunch duration Oj is much larger than the period of the modulation, 
allows one to restrict attention to a small bandwidth around the main mod- 
ulation wavelength at 200 nm. Effectively, we choose a bandwidth between 
167 nm and 251 nm, corresponding to the spectrum f T \ in Fig. 15, which 
is a selection of Fig. 15. Finally, we assume an optics with large numerical 
aperture 0.1 yielding about 10 11 photons. Such number is more than enough 
for the purposes considered in Section 2. 



5 Conclusions 



In this article we presented an extension to our previously proposed method 
for self-seeding H-IPl], which was based on the use of a single Diamond 
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Fig. 15. Squared modulus of the Fourier transform of f x \, |/ T i| 2 in the range between 
167nm and 251nm. 

crystal in Bragg geometry acting as a band-stop filter. Here we extend such 
treatment to a series of two or more crystals tuned at closely-spaced frequen- 
cies, and we demonstrate the possibility of generating doublet spectral lines 
through the self-seeding process. In other words, fully coherent radiation is 
shared between two longitudinal modes. Depending on the crystals settings 
one can have one mode dominating on the other, or a balance between the 
Moreover, the spacing between the two modes can be tuned within 
the FEL gain band, i.e. within 10 eV. Once doublet spectral lines are gen- 
erated, they can be used to study processes dealing with a large change in 
cross-section over a narrow wavelength range. 

Furthermore, possibility of pump-probe experiments with sub-femtosecond 
resolution where an external optical laser is involved are enabled. In fact, a 
power beating between the two closely spaced wavelengths of the doublet 
yields a modulation of the SASE output power on the visible range. By en- 
ergy conservation, such modulation is printed on the electron bunch as an 
energy modulation at the same frequencies, which can be transformed into a 
density modulation with the help of a weak dispersive element. The modu- 
lated electron bunch can then be made radiating coherently with the help of 



2 It should be noted that the idea of radiation shared between longitudinal modes 
should be thought of from a statistical standpoint. In other words, here we are 
talking about an average over ensemble of shots, not about a single realization. 
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any radiator, e.g. an OTR screen. A powerful pulse of coherent radiation in 
the optical range is then produced, which is intrinsically synchronized to the 
X-ray photon pulse, and can be cross-correlated to the external pump-laser. 
Such correlation results in the possibility of measuring the jitter between 
SASE doublet and external optical pump with a few femtosecond accuracy. 
The coherent optical pulse may also be used to monitor the formation of 
the doublet, since it will appear only when the doublet is actually present, 
which enables beating. 

From the technical viewpoint, the present method constitutes a minor up- 
grade of our previously proposed technique, of which it retains flexibility, 
simplicity of implementation, low cost and no risk for the recovering of the 
baseline mode of operation. 
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